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ABSTRACT 


Article history: A solar-powered high temperature differential Stirling engine was considered for optimization using 

Received 28 April 2013 multiple criteria. A thermal model was developed so that the output power and thermal efficiency of 

Accepted 19 June 2013 the so ] ar Stirling system with finite rate of heat transfer, regenerative heat loss, conductive thermal 

bridging loss, finite regeneration process time and imperfect performance of the dish collector could 
be obtained. The output power and overall thermal efficiency were considered for simultaneous maximi¬ 
zation. Multi-objective evolutionary algorithms (MOEAs) based on the NSGA-I1 algorithm were employed 
while the solar absorber temperature and the highest and lowest temperatures of the working fluid were 
considered the decision variables. The Pareto optimal frontier was obtained and a final optimal solution 
was also selected using various decision-making methods including the fuzzy Bellman-Zadeh, LINMAP 
and T0PS1S. It was found that multi-objective optimization could yield results with a relatively low devi¬ 
ation from the ideal solution in comparison to the conventional single objective approach. Furthermore, it 
was shown that, if the weight of thermal efficiency as one of the objective functions is considered to be 
greater than weight of the power objective, lower absorber temperature and low temperature ratio 
should be considered in the design of the Stirling engine. 
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1. Introduction 

Stirling cycle is one of the important primitive standard air cy¬ 
cles of heat engines [1,2]. Suitable efficiency and wide range of 
fuels that can be used for heating can be pointed out as some 
advantages of this type of engine [2-4]. The Stirling engine can the¬ 
oretically be a very efficient engine to convert heat into mechanical 
work at the Carnot efficiency when the ideal regeneration and iso¬ 
thermal compression and expansion processes are considered. 
Thermal limit of the operation of the Stirling engine depends on 
working temperatures on the heater and cooler sides. In most in¬ 
stances, the engine operates with heater and cooler temperatures 
of 923 and 338 K, respectively [5], Engine efficiency ranges from 
30% to 40% coming from a typical temperature range of 923- 
1073 K and the normal operating speed range of 2000-4000 RPM 
[6-11]. Kongtragool and Wongwises [8] studied the influence of 
regenerator efficiency and dead volume on the work as well as effi¬ 
ciency of the machine. However, this study did not include heat 
transfers through the temperature difference at the heat source 
and sink. 
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The idea of coupling for solar concentrators to Stirling engines is 
a new technology which facilitates conversion of solar energy into 
electric power. In this regard, a dish collector with the parabolic 
arrangement of its mirrors is used to concentrate solar radiations 
on the focal point of the collector, in which heat absorber of the en¬ 
gine is located. Therefore, solar energy is collected and concen¬ 
trated towing to a parabola of mirrors. 

In recent years, studying the optimum performance of energy 
systems with thermodynamic models has been an attractive op¬ 
tion for researchers [12-34]. It is believed that the results obtained 
using the FTT lead to a more realistic design of the actual solar en¬ 
ergy systems in comparison to the prior conventional thermody¬ 
namic equilibrium methods. Ladas and Ibrahim [28] defined a 
finite-time parameter as the ratio of working fluid contact time 
to the engine time constant which evolves heat transfer character¬ 
istics of the given design of Stirling engine. They conducted a 
numerical study and plotted variation of the power output versus 
finite-time parameter and change of the power output versus effi¬ 
ciency. Kaushik et al. studied effects of irreversibilities of the 
regenerator and heat transfer process in heat/sink sources [29]. 
Sieniutycz and von Spakovsky [32] generalized thermal exergy to 
the finite-time processes. Sahin and Kodal [33] introduced a new 
thermoeconomic performance analysis based on an objective func¬ 
tion defined as the power output per unit total cost. Tlili investi¬ 
gated effects of regenerating effectiveness and heat capacitance 
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Nomenclature 


A rec absorber (recuperate) area 

A app aperture area 

C concentration ratio 

C v molar specific heat capacity (J mol -1 KT 1 ) 

d deviation index 

d i+ deviation or distance of ith solution from the ideal solu¬ 
tion 

di_ deviation or distance of ith solution from the non-ideal 
solution 

h heat transfer coefficient,(W K 1 or W 1C 4 or W m~ 2 1C 1 ) 

I direct solar flux intensity (W irr 2 ) 

i ith objective 

j jth solution 

n mole number of the working fluid (mol) 

p power (W) 

Q heat transfer (J) 

R universal gas constant (J mor 1 K _1 ) 

T temperature (K) 

W work (J) 

f cyclic period (s) 

k 0 conductive thermal bridge loss coefficient (W K _1 ) 

x temperature ratio of the Stirling engine 


X vector of decision variables 

Greek letter 

X ratio of volume during the regenerative processes (vol¬ 

umetric ratio) 

11 thermal efficiency 

e r effectiveness of the regenerator 

e emissivity factor 

5 Stefan-Boltzman coefficient 

Subscripts 

H absorber (heater) 

HC convection on the high temperature side heat exchanger 
HR radiation on the high temperature side heat exchanger 

L heat sink 

m entire solar-dish Stirling system 

R regenerator 

t Stirling engine 

0 ambient condition, optics 

4-1 process states 


rate of external fluids at the heat source/sink on the maximum 
power and efficiency [34], 

Yaqi et al. developed a mathematical model for the output 
power and overall thermal efficiency of the high temperature dif¬ 
ferential solar dish Stirling engine with the finite heat transfer 
and a non-ideal regenerator and found expressions for the optimal 
absorber temperature which led to the maximized power output of 
the Stirling engine. They evaluated the corresponding thermal effi¬ 
ciency to the maximized power and studied effect of variations in 
engine parameters on the obtained values of the maximized power 
and its corresponding thermal efficiency [17]. They also optimized 
the solar-dish Stirling engine mathematically in a single objective 
optimization approach with only one design parameter (decision 
variable). 

This paper was an extension of the work performed by Yaqi 
et al. [17] with simultaneous consideration of two objective func¬ 
tions instead of the single objective function. In this regard, two 
objective functions including output power and thermal efficiency 
of the entire solar-dish Stirling system were simultaneously con¬ 
sidered. Furthermore, multi-objective optimization was conducted 
with three decision variables including solar absorber temperature 
and the highest and lowest temperatures of the working fluid in 
the heat engine. 

Solving multi-objective optimization problems is a very difficult 
goal which requires simultaneous satisfaction of a number of dif¬ 
ferent and even conflicting objectives. Evolutionary algorithms 
were initially extended and applied during the mid-eighties in an 
attempt to stochastically solve problems of this generic class 
[35]. A reasonable solution to a multi-objective problem is to 
investigate a set of solutions, each of which satisfies the objectives 
at an acceptable level without being dominated by any other solu¬ 
tion [36]. Multi-objective optimization problems generally show a 
possibly uncountable set of solutions, namely Pareto frontier, the 
evaluated vectors of which represent the best possible trade-offs 
in the objective function space. During this term, multi-objective 
optimization of different thermodynamic and energy systems has 
been considered by researchers nowadays [16,37-42]. 

In this paper, a thermodynamic model was developed for 
evaluating the output power and thermal efficiency of a solar-dish 
Stirling engine. The developed model considered convective and 


radiative heat transfers at the hot end of the engine as well as 
the convective heat transfer at the cold end of the engine. More¬ 
over, non-ideal regeneration performance and conductive thermal 
bridge loss between the heat source and heat sink were taken into 
account. The cyclic time of the engine was evaluated based on the 
time spent on the compression, expansion and regeneration pro¬ 
cesses. This model was applied to evaluate the output power and 
thermal efficiency of the entire solar-dish Stirling system for the 
multi-objective optimization process. Considering three design 
parameters (decision variables) including the absorber tempera¬ 
ture and the highest and lowest temperatures of the working fluid 
in the engine, aforementioned objective functions were optimized 
(maximized) using the multi-objective evolutionary algorithm 
(MOEA) and the Pareto frontier was obtained in the objective space 
(power and efficiency space). One class of efficient MOEA, namely 
non-dominated sorting genetic algorithm i.e. NSGA-11, was em¬ 
ployed for maximizing the output power and thermal efficiency 
simultaneously. A final optimal solution from available solutions 
located on the Pareto frontier was selected using three decision¬ 
making approaches including the fuzzy Bellman-Zadeh, LINMAP 
and TOPS1S methods. The final solution obtained by each decision 
maker was compared with the corresponding optimal solutions 
introduced by other decision making tools and the optimal solu¬ 
tion introduced in the single objective optimization as shown in 
[17]. The superiority of solutions obtained using the multi-objec¬ 
tive approach over the corresponding results obtained by the con¬ 
ventional approach ([17]) was discussed and emphasized. 
Sensitivity of the results against variation of the solar-dish Stirling 
engine parameters including the solar collector concentration ra¬ 
tio, effectiveness of the regenerator, conductive thermal bridge loss 
coefficient, volumetric ratio and solar flux intensity was also stud¬ 
ied in detail. 


2. Solar-dish Stirling system 

in solar-dish Stirling systems, a mirror of the parabolic shaped 
concentrator focuses solar light on the focal point of the concentra¬ 
tor where hot end of the Stirling engine is installed. Therefore, solar 
energy with a relatively high temperature is transferred to the hot 
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side heat exchanger of the Stirling engine. Fig. 1 illustrates a sche¬ 
matic for a solar-dish Stirling engine connected to a solar dish con¬ 
centrator. The solar-dish is equipped with a sun tracker which 
tracks the sun in order to have the maximum solar energy transfer 
to the engine when the sun moves during the days. Hence, the 
maximum solar energy is absorbed and transferred to the working 
fluid in the hot space of the engine. 

As shown in Fig. 2, ideal Stirling cycle consists of four processes 
including two isothermal (1-2 and 3-4) and two isochoric pro¬ 
cesses (2-3 and 4-1) in the regenerator. In a real cycle, it is imprac¬ 
tical to have an ideal heat transfer in the regenerator, in which the 
entire amount of absorbed heat (in the process 4-1) is transferred 
to the working fluid in the isochoric heating process (process 2-3). 
Therefore, a heat transfer loss occurs in the regenerator. In addi¬ 
tion, a conductive heat transfer occurs between the heat source 
and sink, namely conductive thermal bridge loss (Qo). 


0_ 




3. Thermodynamic modeling of the solar-dish Stirling engine 


The actual useful heat gain of the dish collector, considering 
conduction, convection and radiation losses, is given as [13]: 

% = IA apP 'h - Arec [h( Th - To) + es(t 4 h - 7t)] (1) 

In Eq. (1 ), / is the direct solar flux intensity, i/ 0 is the collector optical 
efficiency, A rec and A app are the absorber and collector aperture 
areas, respectively, h is the conduction/convection heat transfer 
coefficient in the solar dish, T H and T 0 are the absorber and ambient 
temperatures, respectively, e is the emissivity factor of the collector 
and <5 is the Stefan-Boltzman constant. 

Thermal efficiency i% of the dish collector is obtained as follows 
[17]: 

l s = vti P = llo ~7c [h{T " ~ To) + s5{T " “ T » )] (2) 

It is important to mention that there is a finite heat transfer rate 
in the regenerator, in which the magnitude (Qx) is given by the fol¬ 
lowing expression [17,18,29,30,34]. 

Qr = n C v £ R (T] — T 2 ) ( 3 ) 


where C v is molar specific heat of the working fluid in the regener¬ 
ative processes, e R is effectiveness of the regenerator and n is molar 
mass of the working fluid. Thus, the regenerative heat loss in two 
regenerative processes is given as follows [17,18,29,30,37]: 

AQ r = nCV(l - e*)(T, - T 2 ) (4) 


Owing to the influence of irreversibility of the finite-rate heat 
transfer, time of the regenerative processes is not negligible in 
comparison to that of two isothermal processes [29]. In order to 
calculate time of the regenerative processes, it is assumed that 
temperature of the working fluid in the regenerative processes, 
as a function of time, is given as follows [17-19,30,34]: 



( 5 ) 


where M is proportionality constant which is independent from the 
temperature difference and depends only on the property of regen¬ 
erative material, called regenerative time constant, and ±sign be¬ 
longs to the heating (i = l) and cooling (i = 2) processes, 
respectively [17]. Therefore: 


t3 


T, — T 2 
M, 


( 6 ) 



Total regenerative time is obtained as follows: 

t re = t 3 + U (8) 

For a Stirling cycle, the amounts of heat released to the absorber 
and absorbed by the heat sink are as follows [17]: 

Q.t=lh H c(TH-T,) + h HR (T 4 H -T 4 l )]U 

= nRhLn 1 + nC v ( 1 - £ R )(T, - T 2 ) (9) 

0.2 = h LC (T 2 - T L )t 2 = nRT 2 LnX + nC v ( 1 - s s )(T t - T 2 ) (10) 


where Ch is amount of the heat absorbed by the working fluid at the 
temperature/from the absorber at temperature T H . Similarly, Q 2 is 
amount of the heat released by the working fluid at the tempera¬ 
ture/to the heat sink at temperature T L . In Eq. (9), h HC and h HR are 
convection and radiative heat transfer coefficients at the hot end 
of the engine, respectively. In Eq. (10). h LC is the convective heat 
transfer coefficient on the cold side heat exchanger. Furthermore, 
in Eqs. (9) and (10), R is the universal gas constant and X is 
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volumetric ratio of the engine which is equal to the ratio of volumes 
during the regenerative processes, i.e. 


In addition, there is heat leakage directly from the heat source 
to the heat sink which is directly proportional to the temperature 
difference of sources and the cycle time; i.e.[17,18,29,30], 

Qo = k 0 (T H - T L )t (12) 

where k 0 is conductive thermal bridge loss coefficient, t is cyclic 
period of the temperatures, T H and T L are temperatures of the absor¬ 
ber and heat sink, respectively. 

The net heat absorbed from the heat source and heat released 
Qi to the heat sinkare are obtained using the following expressions: 

Qh = Q,+Qo (13) 

Ql = Q 2 + Qo (14) 

Using Eqs. (6)—(11 ), it is obtained that the cyclic period t is: 


4. Multi-objective optimization 

4.1. Evolutionary algorithms 

Genetic algorithms were developed by John Holland in the 
1960s as a means of importing mechanisms of natural adaptation 
into computer algorithms and numerical optimization [44]. Genet¬ 
ic algorithms apply an iterative, stochastic search strategy to find 
an optimal solution and imitate principles of biological evolution 
in a simplified manner (Fig. 3). In this study, the Pareto frontier 
was obtained using Genetic Algorithm (GA) as a branch of evolu¬ 
tionary algorithm. This method is a powerful optimization tool 
for nonlinear problems [35,37], Pareto optimality is the key con¬ 
cept for establishing a hierarchy among the solutions of a multi¬ 
objective optimization problem in order to determine whether a 
solution is really one of the best possible trade-offs [37]. NSGA-II 
algorithm was utilized here as a multi-objective optimization 
method in order to find the Pareto frontier using genetic algorithm. 
Details of NSGA-II were given in [45]. 


t — tl + t 2 + tre (15) 

_ nRTjLnX + nC„(l - e R )(T, - T 2 ) 
fiHc(T H -Ti) + h HR (T 4 H -Tt) 

nRT 2 LnA-h nCz,(l - £ r )(Ti - T 2 ) f 1 IV 
+ M7W0 + {Mr + W 2 ) {T '- T2) 

Considering cyclic period of the Stirling engine, the power 
output and thermal efficiency of the engine are given as 


follows: 


p t t 

(16) 

n =^h-& 
nt du 

(17) 


Substituting Eq. (15) in Eqs. (16) and (17) leads to the fol¬ 
lowing expressions for power and thermal efficiency of the Stir¬ 
ling engine: 


Tj - T 2 


Ti+M(T|-r 2 

cUh 


+M(T,-T 2 ) 7 2 +M(T,-T 2 ) . c (T T \ 

b)+MrS-rj) + hlc(T2 - Tl) 


(18) 


h,= 


t,-t 2 


r 1+ M(r 1 -r 2 ) + [k 0 (r H -T t )][j 


Ti+MgirS) , t 2 +M(T,-t 2 ] 
hHc(r H -T,)+h H ,er 4 H -T‘) '■icin 


i^ + F,(T,-T 2 )j 
(19) 


where 


M = 

F, = 


C„( 1 - £r) 

RLnX 

1 / 1 1 \ 

nRLnX \Nh + U 2 ) 


(20a) 

(20b) 


Thermal efficiency of the system is a product of thermal effi¬ 
ciency of the collector (Eq. (2)) and optimal thermal efficiency of 
the Stirling engine (Eq. (19)) [17]. Therefore: 


'lm = >h>h 


>1m = {>1o-l{h(TH-T 0 ) + es(Tt 1 -Tt)}} 


(21a) 


(21b) 


l+M(l-x)+[k 0 (T H -r L )] 


{ h HC^ T H~ T \opt> +l, HR 


1+M(l-x) _ [ x+M(l-x) 

h Lc( xT 1opt~ T 0 


( i 4_ t 4 ) 

^ H lopt J 


(1 +b) 


4.2. Objective functions, decision variables and constraints 

Two objective functions of this study included output power of 
the system (P) and the system thermal efficiency (f/ m ) denoted by 
Eqs. (18) and (21b), respectively. 

In this paper, three decision variables were considered as follows: 

T,: highest temperature of the working fluid at isothermal pro¬ 
cess (process 3-4 in Fig. 1). 

T 2 : lowest temperature of the working fluid at isothermal pro¬ 
cess (process 1-2 in Fig. 1). 

T h : absorber temperature. 

It is required to mention that the number of decision variables 
could be increased to four decision variables if volumetric ratio of 
the engine (7.) is considered an additional decision variable. How¬ 
ever, in this research, it was required to have the highest possible 
value of X in order to have the highest values for the maximum val¬ 
ues of the output power and thermal efficiency. In other words, X is 
always selected at its upper bound values defined for a variation of 
this variable. Due to the engineering limitation, X is usually chosen 
from a range between 2.0 and 5.0. Therefore, if this constraint is 
implemented for X, optimization leads to X = 5.0 for this variable. 
Thus, the optimization would be meaningless for this variable 
and it is selected at its highest possible value. Hence, in this paper, 
it was removed from the list of decision variables. 

The objective functions were solved with respect to the follow- 


ing constraints: 


Tj > T 2 

(22) 

T 2 > 320 °C 

(23) 

700 °C^T h ^ 1600 °C 

(24) 

0.4 < T 2 /T, ^ 0.7 

(25) 



Fig. 3. Schematic of the objective space. 
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5. Decision-making in multi-objective optimization 

In multi-objective optimization, a process of decision-making is 
required for selection of the final optimal solution from available 
solutions. There are several methods for decision-making process 
in decision problems. These methods can be employed for selection 
of a final optimal solution from the Pareto frontier. Since dimension 
of various objectives in a multi-objective optimization problem 
might be different (for example, in the present case, thermal effi¬ 
ciency objective had no dimension while dimension of the output 
power was in Watt), before any decisions, dimensions and scales of 
objective space should be unified. In this regard, objective vectors 
should be non-dimensioned before the decision-making process. 
There are several methods of non-dimensioning which are utilized 
in decision making and include linear non-dimensioning, Euclidean 
non-dimensioning and fuzzy non-dimensioning. In this fuzzy Bell¬ 
man—Zadeh, LINMAP and TOPSIS methods of decision making meth¬ 
ods were employed and the Bellman-Zadeh method utilized fuzzy 
non-dimensioning while other methods (LINMAP and TOPSIS) em¬ 
ployed Euclidean non-dimensioning. Details of fuzzy non-dimension¬ 
ing and Euclidean non-dimensioning are given in [43], 

5.1. Bellman-Zadeh decision-making method 

A final decision is defined by the Bellman and Zadeh model as 
the intersection of all fuzzy criteria and constraints and is repre¬ 
sented by its membership function. In this method, a matrix of 
membership function is defined. The column of this matrix in¬ 
cludes fuzzy membership function for each objective. Rows of 
the membership matrix reflect magnitudes of the membership 
functions for each solution of the Pareto frontier. Therefore, the 
number of columns in the membership function matrix is equal 
to the number of objectives and the number of rows is equal to 
the number of solutions located on the Pareto frontier. Details of 
the method of defining the membership function were presented 
by Sayyaadi and Mehrabipour [43]. The fuzzy membership func¬ 
tion for each solution is set at the minimum value of the member¬ 
ship functions of all objectives for the proposed solution. Hence, a 
vector of membership is obtained, in which the elements are the 
minimum membership function of objectives at each solution. Fi¬ 
nally, maximum of these minimums, i.e. an element of the mem¬ 
bership vector with maximum values of membership function, is 
selected as a final solution. In other words, in the Bellman-Zadeh 
fuzzy decision making, a maximum of minimums (minimum 
membership function of objectives) is selected as the final optimal 
solution. Readers are encouraged to refer to [43] for more details. 

5.2. LINMAP decision-making method 

An ideal point on the Pareto frontier is a point at which each 
objective is optimized regardless of the satisfaction of other objec¬ 
tives. It is clear that, in the multi-objective optimization, it is 
impossible to have each objective in its optimal condition which 
can be obtained in a single objective optimization. Therefore, the 
ideal point is not located on the Pareto frontier. In the LINMAP 
method, after Euclidean non-dimensionalization of all objectives, 
the spatial deviation of each solution on the Pareto frontier from 
the ideal point is determined. Finally, a solution with a minimum 
spatial distance from the ideal point is selected as a desired final 
optimal solution. More details of the LINMAP method were pre¬ 
sented by Sayyaadi and Mehrabipour [43]. 

5.3. TOPSIS decision-making method 

In this method, beside the ideal point, a “non-ideal point” is de¬ 
fined. The non-ideal point is the ordinate in objective space, in 


which each objective has its worst value. Therefore, beside the 
solution spatial distance from the ideal point, the solution spatial 
distance from the non-ideal one is used as a criterion for selecting 
the final solution. In the TOPSIS method, a solution with the max¬ 
imum distance from the non-ideal point and minimum distance 
from the ideal point is selected as a desired final optimal solution. 
More details of the TOPSIS method could be found in [43]. 

6. Result and discussion 

Output power and thermal efficiency of the solar-dish Stirling 
engine were maximized simultaneously using the multi-objective 
optimizing method, which worked based on the NSGA-II algo¬ 
rithm. In this regard, optimization was performed by objective 
functions, as expressed by Eqs. (18) and (21b), and constraints, as 
stated in Eqs. (22)-(25). Design parameters (decision variables) 
of the optimization were temperature of the working fluid in the 
high temperature isothermal process (Ti), temperature of the 
working fluid in the low temperature isothermal process (T 2 ) and 
the absorber temperature (T H ). 

To be consistent with the previous works, specifications of the 
solar-dish Stirling engine were considered as follows [17]: 

x = 0.5, h HC = h LC = 200 W K 1 , h HR = 4 x 10~ 8 W I< 4 , n 

= 1 mol, 2 = 2, s R = 0.90, e = 0.90, R = 4.3 J moP 1 K \C„ 

= 15 J mor 1 K 1 , <5 = 5.67 x lO^Wm- 2 I< 4 ,T t = 320 K, T 0 

= 300 K, h = 20\N nr 2 K 1 . T H = 1100 K, / = 1000 W nr 2 , C 

= 1300. (^ + = 2 x 10- 5 s K 1 

Vjw, m 2 J 

Fig. 4 illustrates the Pareto frontier in the proposed objectives’ 
space obtained by the multi-objective evolutionary algorithm 
(NSGA-II). In this figure, the ideal and non-ideal are specified. Fur¬ 
thermore, three final solutions selected by the fuzzy, LINMAP and 
TOPSIS decision makers are indicated in Fig. 4. 

Table 1 compares optimal results obtained in this paper using 
three decision making methods with the corresponding results ob¬ 
tained in a single objective-single decision variable work [13], Fur¬ 
thermore, this table compares solutions with the ideal solution 
(ideal solution is illustrated in Fig. 4). It is clear that thermal effi¬ 
ciency of the final solutions selected using the TOPSIS, LINMAP 
and fuzzy decision-makers was 4.87%, 4.13% and 1.23% less than 
the corresponding thermal efficiency obtained in [17], respectively. 
Thermal efficiency of the solutions obtained using the TOPSIS, LIN¬ 
MAP and fuzzy decision makers was 8.06%, 7.32% and 4.42% less 



8 10 12 14 16 18 20 22 24 


P(kW) 

Fig. 4. Pareto optimal frontier in the objective space. 
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Table 1 

Comparison of the final optimal solutions of this paper (multi-objective, multi-variable) with the corresponding results obtained in [17] (single objective-single variable 
optimization). 


Optimal solution 

Decision making method 

Decision variable 

r, (K) r 2 (i<) 

Th 00 

Objectives 

P(W) 

dm 

Deviation index from the ideal solution 

d 

This work 

TOPSIS 

1248.3 

571.4 

1565.6 

22286.8 

0.2594 

0.0682 


LINMAP 

1248.0 

539.4 

1569.2 

21587.4 

0.2668 

0.0650 


FUZZY 

1192.7 

506.8 

1449.5 

18113.8 

0.2958 

0.0786 

Ref. 113] 

- 

923.9 

462.0 a 

1100.0 3 

10164.2 

0.3081 

0.2252 

Ideal solution 

- 

- 

- 

- 

22800.0 

0.3400 

0.0000 


a This parameter was not considered a decision variable in Ref. [17]. Its fix value in [17] is indicated in the table. 



4 6 8 10 12 14 16 18 20 

P (kWatt) 


Table 2a 


Optimal solutions of the Stirling engine at 
(£r = 0.90). 

various 

volumetric ratios of the engine 

X 

Decision making 
method 

Decision variable 

Objectives 

h (i<) 

r 2 (K) 

Th (I<) 

P(W) 

dm 

2.0 

TOPSIS 

1248.3 

571.4 

1565.6 

22286.8 

0.2594 


LINMAP 

1248.0 

539.4 

1569.2 

21587.4 

0.2668 


FUZZY 

1192.7 

506.8 

1449.5 

18113.8 

0.2958 

3.0 

TOPSIS 

1250.5 

539.8 

1589.9 

26356.8 

0.2865 


LINMAP 

1250.5 

539.8 

1589.9 

26356.8 

0.2865 


FUZZY 

1208.4 

492.6 

1465.2 

21410.8 

0.3279 

4.0 

TOPSIS 

1275.1 

544.9 

1593.0 

28726.0 

0.2976 


LINMAP 

1273.3 

538.4 

1570.4 

27832.5 

0.3062 


FUZZY 

1226.7 

509.9 

1436.7 

22312.1 

0.3427 

5.0 

TOPSIS 

1261.3 

534.5 

1575.2 

29002.6 

0.3102 


LINMAP 

1261.6 

534.7 

1548.9 

28133.4 

0.3178 


FUZZY 

1216.8 

498.2 

1443.1 

23254.3 

0.3509 



Fig. 5. Pareto optimal frontier in objective space at various (a) degrees of 
effectiveness of regenerator; (b) volumetric ratio of the engine. 


than that of the ideal solution, respectively. This table indicates 
that the best value for the thermal efficiency obtained in this paper 
was 1.23% less than that of [17] and 4.42% less than that of the 
ideal solution. Comparison of the output power of the engine ob¬ 
tained using three decision making with the corresponding output 
power obtained in [17] indicated that output power of the solu¬ 
tions obtained by TOPSIS, L1NMAP and fuzzy methods was 
119.3%, 112.4% and 78.2% higher than the output power obtained 
in [17]. The output power obtained in this paper for three solutions 
using the TOPSIS, LINMAP and fuzzy decision making methods was 
2.2%, 5.3% and 20.6% less than the output power of the ideal solu¬ 
tion, respectively. Results in Table 1 indicate that, although ther¬ 
mal efficiency of solutions obtained in this work was less than 
the corresponding thermal efficiency obtained in the single objec¬ 


tive work [17], much greater enhancement was obtained for the 
output power objective in comparison to [17]. In other words how¬ 
ever, the least thermal efficiency of the Stirling engine obtained in 
this work (using TOPSIS decision maker) was 4.87% less than the 
corresponding thermal efficiency of the single objective-single var¬ 
iable work [17] and the least output power (obtained by the fuzzy 
decision maker) was 78.2% higher than the corresponding output 
power of [17]. For better insight into the status of various solutions 
obtained here in comparison to the solution of the single objective- 
single variable work [17], a deviation index characterizing devia¬ 
tion of each solution from the ideal solution was also introduced. 


d + = \J(Pn ~ Pn,ideal? + (»7m„ “ Vm^deai? ■ 

d - — yj(Pn ~ Pn,non-ideal) “F Im n ~ *1 m n ,non-ideal) 


d = 


(d+) +( d -) 


(26) 

(27) 

(28) 


p n and /7 m „ denote Euclidian non-dimensioned output power and 
thermal efficiency (indexes ideal and non-ideal denote the corre¬ 
sponding values of objective functions at the ideal and non-ideal 
points, respectively). The last column of Table 1 denotes deviation 
index (d) for the solution of each optimization approach. As is clear, 
deviation indexes for the multi-objective optimization approach of 
this paper were 0.0682, 0.0650 and 0.0.0786 for the solutions ob¬ 
tained by TOPSIS, LINMAP and fuzzy decision making methods, 
respectively. These values were comparable with the corresponding 
value obtained for the single objective-single variable optimization 
[17], where deviation index is 0.2252. This issue implies that devi¬ 
ation indexes of the multi-objective, multi-variable optimization for 
all decision making methods were much lower than the single 
objective-single variable optimization (the ideal deviation index 
was zero). Also, the last column of Table 1 indicates that the LIN- 
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Table 2b 

Optimal solutions of the Stirling engine at various degrees of effectiveness of the 
regenerator (2 = 2.0) parameter. 


e R 

Decision making 
method 

Decision variables 

Objectives 

T, (K) 

T 2 (K) 

M K) 

P(W) 

r\m 

0.5 

TOPSIS 

1195.5 

692.2 

1535.6 

12132.6 

0.1289 


LINMAP 

1199.7 

634.4 

1539.6 

11536.1 

0.1348 


FUZZY 

1154.1 

582.6 

1411.0 

9376.0 

0.1507 

0.6 

TOPSIS 

1222.9 

655.8 

1575.6 

13845.9 

0.1453 


LINMAP 

1218.2 

651.3 

1557.6 

13556.8 

0.1480 


FUZZY 

1169.4 

575.3 

1420.4 

10649.2 

0.1706 

0.7 

TOPSIS 

1240.8 

674.7 

1555.7 

15853.0 

0.1671 


LINMAP 

1220.7 

628.9 

1543.2 

15153.9 

0.1745 


FUZZY 

1167.8 

539.4 

1453.6 

12262.3 

0.1965 

0.8 

TOPSIS 

1243.3 

601.2 

1577.3 

18281.2 

0.2047 


LINMAP 

1236.7 

585.0 

1529.4 

17112.3 

0.2164 


FUZZY 

1189.6 

527.6 

1446.1 

14422.1 

0.2372 

0.9 

TOPSIS 

1248.3 

571.4 

1565.6 

22286.8 

0.2594 


LINMAP 

1248.0 

539.4 

1569.2 

21587.4 

0.2668 


FUZZY 

1192.7 

506.8 

1449.5 

18113.8 

0.2958 


MAP decision making method led to a lower value for deviation in¬ 
dex; hence, in this paper, the solution that was selected using the 
L1NMAP decision making was considered a final desired optimal 
solution of the multi-objective, multi-variable optimization for the 
solar-dish Stirling engine. It should be noted that this is not a gen¬ 
eral rule; as in other cases, other decision making methods may 
yield better results. 



T h (K) 



t h (K) 


Fig. 5a and b shows the Pareto optimal frontier at various de¬ 
grees of effectiveness of the regenerator and volumetric ratios of 
the engine, respectively. 

Optimal solutions introduced by three decision makers in the 
cases that the engine was designed with different volumetric ratios 
are indicated in Table 2a. Similar results at various degrees of effec¬ 
tiveness of the regenerator are indicated in Table 2b. As expected, 
engines with higher values of volumetric ratio and regenerator 
effectiveness had higher optimal values for both objectives (power 
and thermal efficiency). 

Fig. 6a-d illustrate absorber temperature of the Stirling engine 
for solutions of the Pareto frontier at various degrees of effective¬ 
ness of the regenerator (e R ) and volumetric ratio (2). Fig. 6a and b 
shows that, for solutions of the Pareto frontier with higher values 
of thermal efficiency, the absorber temperature (T H ) was lower. 
As previously mentioned, all the solutions located on the Pareto 
frontier could be considered a potential final optimal solution of 
the multi-objective optimization. In this regard, the solutions lo¬ 
cated on the left side of the Pareto frontier (Fig. 4) were solutions 
with higher thermal efficiency and lower output power. It means 
that, if a weight factor were considered for each objective (a value 
of ^1.0), solutions located on the right side of the Pareto frontier 
would have higher weight factors for the thermal efficiency and 
lower weight factors for the output power. Fig. 6a and b imply that, 
for solutions of the Pareto frontier with a higher weight factor for 
thermal efficiency, values of absorber temperature were lower. 
When the weight factor for thermal efficiency was 1.0 (the first 
solution on the left side of the Pareto frontier), the absorber tem¬ 
perature was 1100K, which was consistent with the results ob- 



T „(K) 



T h (K) 


Fig. 6. Optimal solutions of the Pareto frontier as a function of the absorber temperature (7 H ), effectiveness of the regenerator (£ R ) and volumetric ratio of the engine ( 2 ) (a) 
dm - Th - Or, (b) ;/ m - T h - 2 ; (c) P - T„ - e R ; (d) P - T H - 2 . 
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Fig. 7. Optimal solutions of the Pareto as a function of the temperature ratio of the engine (x), effectiveness of the regenerator (e R ) and volumetric ratio of the engine (7) (a) 

rim-x- e R ; (b) t; m - x - 2; (c) P - x - e B ; (d) P - x - X. 


Table 3 

Sensitivity of objective functions and decision variables to variations of the system parameter. 


Parameters 

Proposed value 

Possible variations 

Power (kW) 

Efficiency 

T, (K) 


T 2 (I<) 


Th (I<) 


Lower limit 

Upper limit 

a 

b 

a 

b 

a 

b 

a 

b 

a 

b 

X 

2 

2 

5 

17.83 

23.09 

0.30 

0.35 

1225.1 

1212.9 

500.1 

486.1 

1454.7 

1449.2 

£ r 

0.9 

0.5 

0.9 

9.24 

17.25 

0.15 

0.30 

1170.3 

1196.2 

558.7 

515.1 

1437.4 

1410.2 

7 (W irr 2 ) 

1000 

500 

1500 

9.24 

18.52 

0.15 

0.32 

1203.9 

1216.2 

567.4 

501.5 

1433.3 

1479.0 

c 

1300 

700 

2000 

16.63 

19.18 

0.25 

0.32 

1134.6 

1231.8 

513.9 

514.7 

1370.7 

1491.8 

M WK->) 

2.5 

1.5 

5 

17.52 

18.72 

0.31 

0.28 

1220.5 

1227.1 

503.1 

506.6 

1438.5 

1482.2 


a Value corresponds to the lower limit of the proposed parameter. 
b Value corresponds to the upper limit of the proposed parameter. 


Table 4 

Error analysis based on the mean absolute percent error (MAPE) method. 


Error analysis 

FUZZY 


LINMAP 


TOPSIS 



P 

y] m 

P 

r\m 

P 

V]m 

Maximum error (%) 

2.02 

3.65 

1.71 

4.37 

2.70 

3.82 

Average error (%) 

0.64 

1.65 

1.14 

1.90 

1.33 

1.54 


tained in [13]. The absorber temperature increased to 1600 I< when 
the weight factor for thermal efficiency decreased to 0.0; i.e. 
weight factor for the output power was considered 1.0 (the first 
solution on the right side of the Pareto frontier). Therefore, the ab¬ 
sorber temperature of a final solution with greater weight for the 
output power should be kept at a higher value than the corre¬ 


sponding absorber temperature of solutions with higher weight 
factor for the thermal efficiency. This fact is also elucidated in 
Fig. 6c and d, which present the output power of the Pareto solu¬ 
tions as a function of the absorber temperature. In addition, 
Fig. 6a-d demonstrate that output power and thermal efficiency 
increase with the incraese in effectiveness of the regenerator and 
volumetric ratio of the engine. 

Fig. 7a-d represent similar graphs for solutions of the Pareto 
frontier. These figures illustrate the magnitude of objectives for 
solutions of the Pareto frontier as functions of the temperature ra¬ 
tio (x), effectiveness of the regenerator (e R ) and volumetric ratio of 
the engine (a). They also show that, for solutions with higher val¬ 
ues of weight factors for thermal efficiency, temperature ratio of 
the engine was lower than the corresponding temperature ratio 
of the solutions with lower weight factors for thermal efficiency. 
Flence, when a higher weight factor was considered for the output 
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power, the engine should be designed with a higher temperature 
ratio. 

6.3. Sensitivity analysis 

Optimization is usually performed while some parameters are 
considered to be fixed. In sensitivity analysis, effect of variations 
in these kinds of parameters on the final selected solutions is stud¬ 
ied. Sensitivity analysis was performed here in order to evaluate ef¬ 
fect of variations in system parameters including volumetric ratio, 
effectiveness of the regenerator, solar flux intensity, concentration 
ratio of the solar dish and conductive thermal bridge loss coeffi¬ 
cient. Table 3 indicates effect of variations in aforementioned 
parameters of objective functions and decision variables. 

This table implies high sensitivity of objective functions against 
variations of the proposed system parameters. Therefore, in evalu¬ 
ation of the parameters including volumetric ratio, effectiveness of 
the regenerator, solar flux intensity, concentration ratio of the solar 
dish and conductive thermal bridge loss coefficient, special care 
must be paid and results obtained in the condition in which these 
parameters are evaluated should be generalized to other cases 
with caution. 

6.2. Error analysis 

In this paper, error analysis was performed based on the mean 
absolute percent error (MAPE) method. In this regard, 30 runs for 
the computer code were performed and a final solution for each 
run was obtained using the fuzzy, LINMAP and TOPSIS decision 
making methods. Then, magnitudes of each objective (output 
power and thermal efficiency) were compared with its correspond¬ 
ing best value (presented as a final solution in this paper). The first 
row of this table represents average error of each objective for 
solutions of three decision making methods. The last row of this ta¬ 
ble indicates a maximum error of each objective obtained using 
three decision making methods. 

According to the data in Table 4, it can be found that average er¬ 
ror of the solutions obtained by three decision making methods 
were 1.04% and 1.70% of the output power and thermal efficiency, 
respectively. Similarly, it can be found that the maximum error of 
solutions obtained by three decision making methods were 4.37% 
and 2.02% of the output power and thermal efficiency, respectively. 

7. Conclusions 

A thermodynamic model was applied to determine output 
power and thermal efficiency of the solar-dish Stirling heat engine. 
The output power and thermal efficiency of the engine were con¬ 
sidered simultaneously for multi-objective optimization, in which 
temperature of the working fluid in the high temperature isother¬ 
mal process, temperature of the working fluid in the lower temper¬ 
ature isothermal process and the absorber temperature were 
considered design parameters. A final optimal solution was se¬ 
lected from available solutions of the Pareto frontier using three 
decision making methods including LINMAP, TOPSIS and fuzzy 
methods. It was shown that the multi-objective, multi-variable ap¬ 
proach could lead to a more desirable design of the Stirling cycle if 
compared with the single objective-single variable approach per¬ 
formed in the previous researches. The priority of the multi-objec¬ 
tive, multi-variable approach over the traditional single objective- 
single variable one was proven in all cases and significant lower 
deviation indexes were obtained for the multi-objective, multi- 
variable approach. Moreover, in this case, it was observed that 
the LINMAP decision making method could yield a final solution 
with the least deviating index from the ideal solution. The volu¬ 


metric ratio effect and effectiveness of the regenerator on the final 
selected solution was studied and it was found that a final optimal 
solution with higher values of objectives (thermal efficiency and 
output power) could be achieved if these parameters were in¬ 
creased. It was also shown that, in the multi-objective optimization 
of the solar-dish Stirling engine, if a great weight factor is consid¬ 
ered for the thermal efficiency, absorber temperature and temper¬ 
ature ratio of the engine should be considered at a lower value in 
comparison to the case in which a lower weight factor is consid¬ 
ered for the thermal efficiency. Relatively high sensitivity of the fi¬ 
nal selected optimal solution against variation of the system 
parameters including the volumetric ratio, effectiveness of the 
regenerator, solar flux intensity, concentration ratio of the solar 
dish and conductive thermal bridge loss coefficient was found; 
therefore, evaluation of these parameters in the design process of 
the solar-dish Stirling should be performed precisely. 
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